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We investigate whether there are any cosmological evidences for a scalar field with a mass and cou- 
pling to matter which change accordingly to the properties of the astrophysical system it "lives in" , 
without directly focusing on the underlying mechanism that drives the scalar field scale-dependent- 
properties. We assume a Yukawa type of coupling between the field and matter and also that the 
scalar field mass grows with density, in order to overcome all gravity constraints within the solar 
system. We analyse three different gravitational systems assumed as "cosmological indicators" : su- 
pernovae type la, low surface brightness spiral galaxies and clusters of galaxies. Results show that: 
a. a quite good fit to the rotation curves of low surface brightness galaxies only using visible stellar 
and gas mass components is obtained; b. a scalar field can fairly well reproduce the matter profile 
in clusters of galaxies, estimated by X-ray observations and without the need of any additional 
dark matter; c. there is an intrinsic difficulty in extracting information about the possibility of a 
scale-dependent massive scalar field (or more generally about a varying gravitational constant) from 
supernovae type la. 
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I. INTRODUCTION 

Present cosmological and astrophysical observations 
clearly depict a universe dominated by "dark" compo- 
nents, being it made by dark matter for 25% and by 
dark energy for 71%, while ordinary baryonic matter con- 
tributes for only the remaining 4% [s^ . 

Dark matter has a long history, being "introduced" by 
[9^ : for solving the problems of high mass-to- light ra- 
tios of galaxy clusters and of the rotation curves of spiral 
galaxies. Later on, several versions of the the so-called 
cold dark matter model (CDM) have been built starting 
from the assumption that a large amount of non-baryonic 
matter, matter non-interacting with the electromagnetic 
radiation but only detectable by its gravitational inter- 
action with visible matter, could account for the obser- 
vations in the framework of the standard Newtonian dy- 
namics. But even if its clustering and distribution prop- 
erties are fairly well known at every scale (see [56j for the 
most used model), its nature is unknown, up to now, at 
a fundamental level. 

Dark energy has a more recent history: it was "in- 
troduced" about ten years ago while reconstructing the 
Hubble diagram of Type la Supernovae (SNela) obser- 
vations, from which it was deduced that the universe 
is now accelerating [l], [2^ [g^, [zll- While growing the 
quantity of available cosmological data (measurements 
of cluster properties as the mass, the correlation func- 



tion and the evolution with redshift of their abundance 
H, 0, m, HI]; the already mentioned Hubble diagram 
of SNela; the optical surveys of large scale structure 
[27I [ssl . Issj : the anisotropies in the cosmic microwave 
background (29l . [sol [78| : the cosmic shear measured from 
weak lensing surveys |69l. [STj and the Lyman - a forest ab- 
sorption |28l [53} ) ■ more evidences towards a spatially flat 
universe with a subcritical matter content and undergo- 
ing a phase of accelerated expansion have been collected. 
Despite of the existence of universe acceleration has 
been clearly undisclosed, the nature and the fundamental 
properties of the underneath physical mechanism remain 
essentially unknown notwithstanding the great theoreti- 
cal efforts made up to now. By simply adding a constant 
to the dynamical equations of universe, the cosmological 
constant A [l^ [t^ , in the context of CDM model, it was 
defined a new model, the ACDM, which quickly became 
the consensus model because it provides a good fit to 
most of the data [tI, [t^, [13 giving a reliable snapshot of 
the today observed universe. Nevertheless, it is affected 
by serious theoretical shortcomings that have motivated 
the search for more general and alternative candidates 
generically referred to as dark energy. Such models range 
from scalar fields rolling down self interaction potentials 
to phantom fields, from phenomenological unified mod- 
els of dark energy and dark matter to alternative gravity 
theories [13, M, MM,^- 

In the last three decades, scalar fields have played an 
important role in both cosmology and particle physics 
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[5|, l52| . Scalar fields have been postulated as means to 
explain the early and late time acceleration of the uni- 
verse. However, it is almost always the case that such 
fields interact with standard matter: either due to a di- 
rect Lagrangian coupling or indirectly through a coupling 
to the Ricci Scalar or as the result of quantum loop cor- 
rections. Both for inflation in the early universe and for 
dark energy, such couplings can lead to problems. In in- 
fiation, for example, couplings might destroy the flatness 
of the potential needed to drive a period of inflation. If 
there are scalar fields which permeate the universe to- 
day and have non-zero couplings to matter, then they 
would induce an additional force in nature. If the scalar 
field self-interactions are negligible, then the experimen- 
tal bounds on such a field are very strong: either the 
couplings to matter are much smaller than gravity, or 
the scalar fields arc very heavy, so that they mediate a 
short-ranged interaction. 

A certain class of theories have been proposed, in 
which the scalar field properties depend on the environ- 
ment: its mass depends on the local environmental den- 
sity. These are the so called Chameleon Field Theories, 
proposed by [i^, that employs a combination of self- 
interaction and couplings to matter of the scalar field to 
avoid the most restrictive of the current gravity bounds. 
In these models a scalar field couples to matter with 
gravitational strength, in harmony with general expec- 
tations from string theory, whilst, at the same time, re- 
maining relatively light on cosmological scales. It was 
found that local gravity constraints are (roughly) satis- 
fied as long as the mass-scale of the potential satisfies 
M < (lmm)~^. This coincides with the scale associ- 
ated with the late time acceleration of the universe, and 
it is surprising that it should come from local experi- 
ments. Chameleon models have been subject to many 
studies, from laboratory exp eriments up to cosmologi- 
cal probes ^M, M,\Mm,\M, MM, tH-il ■ Lately, 
chameleons models where also studied in the context of 
Modified Gravity, in particular the so called /(7?)-gravity 
pll [l9l . Isoj . In such a case, good results of general rela- 
tivity, according to the standard probes, are reproduced 
at local scales while dark energy (accelerating behavior) 
and dark matter (dynamical clustering) effects are repro- 
duced at larger scales. The major issue is to select suit- 
able /(i?)-models capable of matching the density pro- 
files at the various gravitational scales, as discussed e.g. 
in [3. 

Chameleon fields evade tight gravity constraints via 
the so called chameleon mechanism: it consists on a field- 
generated force, a sort of fifth force, to become short- 
ranged in highly dense regions, and long-ranged in low 
density regions. Such feature would imply that at dif- 
ferent astrophysical scales the fifth force felt by matter 
would be suppressed or enhanced according to the local 
astrophysical density. 

In this work we investigate whether there is evidence 
for a coupling between baryonic matter and a massive 
scalar field which could mimic and replace the con- 



tributions of a possible dark matter component. In 
particular, motivated by chameleon models, we aim 
to verify if it is possible to observationally detect a 
scalar field whose mass (or, equivalently, interaction 
length) and coupling may change with scale, matching 
different astrophysical observations. In this case, by 
astrophysical observations, we mean data from SNela, 
low surface brightness (LSB) dwarf galaxies and, fi- 
nally, clusters of galaxies. The range of scales is very 
wide and we are adopting photometric and spectro- 
scopic data to probe the mechanism at different redshifts. 

We want to stress here that we are just interested to 
investigate whether the data would favor a model where 
baryons may be coupled to a scalar field whose mass and 
coupling change with scale. If such a change is due to a 
chameleon mechanism (associated to a local-density vari- 
ation) or due to some other mechanism, that is not the 
main concern in this article and so we will not compute 
specific predictions for any particular model. 

Chameleon models arc highly non linear so that, in 
principle, no superposition principle for such a non-linear 
field would be possible. Due to that, it is extremely 
difficult to build the gravitational potential of an ex- 
tended astrophysical system as the ones we analyse in 
this work. Fortunately, we are not specifically investi- 
gating chameleon models; we are here studying a general 
coupled scalar field model and asking the data what are 
the preferred values of its mass and coupling at different 
scales. In order to compute an extended gravitational 
potential (in a non-linear case as the chameleon model 
would be), we would need to properly study the non- 
linear regime of structure formation within this type of 
models. But this is well beyond the scope of this article. 
In spite of being of utmost importance to set up N-body 
simulations with heavy scalar field models, such task is 
extremely complex. In fact up to nowadays most of the 
N-body simulations assume light scalar fields which do 
not cluster at small scales. That is obviously not the case 
of chameleon fields and such investigation is not consid- 
ered in the present paper. 

The article is organized as follows: In § (|n| we give a 
brief but exhaustive summary of all the main properties 
of the scalar field theory. In § piip we accurately de- 
scribe the used astrophysical and cosmological data and 
the theoretical model defined for any of them. In § (|IV[) 
we expose our results with a discussion on their impli- 
cations for a more general and comprehensive theory of 
gravity. Conclusions are drawn in § (|V|) 



II. THE SCALAR FIELD THEORY 

A general action governing the dynamics of the 
chameleon (scalar) field <f> can be of the form: 

|d^xy^|^7^-i(a<^)^-l^(0)| 
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Figure 1. Chameleon potential in low and high local density environment. 



- I d4x£™(^«,5;.), (1) 

where g is the determinant of the metric g^j,. TZ is the 

Ricci scalar, ^pm are the various matter fields and £,„ is 

the Lagrangian density of ordinary matter. 

In the expression for the reduced Planck Mass, Mpi = 

(SttG*)"^/^, is the bare gravitational constant and 

differs from the usually measured one. 

This can be better understood if we consider the most 

general case of a scalar-tensor model, whose action (Eq. [T] 

is a particular case of this) is [37l |: 

-Vm-Jd<xC,n{^l^\g^^,). (2) 

/(i?)-gravity models too can be enclosed in this general 
case since it is straightforward to show that, by suitable 
manipulations, they are a subclass of Eq. ([2]); see for 
example [2^. From this action one can naively define 
Newton's gravitational constant as: 

G 

Gw = — , (3) 

but Gjv has not the same physical meaning as Newton's 
gravitational constant in general relativity. The actual 
Newtonian force measured in a Cavendish-type experi- 
ment between two test masses will experience an effective 
coupling constant 

G,ff = ^{l + am = GN{l + am. (4) 

where the term G^,/F is due to the (average) exchange 
of gravitons between the two bodies, while G^/F ■ a{(j)) 
comes from the exchange of a scalar particle between 
them, with the analytical expression for a depending on 
the particular scalar theory one considers. It is also clear 
from this expression, that what in general relativity is a 
true constant, now becomes a possible function of time 



and radius, so the use of the term "constant'' is quite in- 
appropriate. When we will speak about a varying grav- 
itational constant we always refer to the effective scalar 
field gravitational constant expression. 

The key ingredient in studying the cosmological dy- 
namics of a field, when it supports a chameleon mecha- 
nism, is that the scalar field feels both a potential, V{(j)), 
and a coupling to matter depending on p, the local den- 
sity of matter, and on the coupling constant, /3. At the 
end the field dynamics is governed by an effective poten- 
tial, 

y,/^(</.) = l/(0)+^p,e^'<^/^- . (5) 

i 

If V{(t)) is a run-away potential and /3i > 0, the effective 
potential has a minimum at (f>min satisfying the condi- 
tion: 

K0(0™n) + E 1^P^<^^'*^'''' = , (6) 

Mpi 

and the effective mass of the field of small perturbations 
about the minimum is 

= = VM^„..n) + E ■^P^e^^''/''-^ . (7) 

So the effective potential depends on the matter density 
and both the minimum in potential and the mass of the 
scalar field are function of the local density, as shown in 
Fig. ([1]). As density increases, the minimum in potential 
shifts to smaller values of the field and the mass of small 
fiuctuations increases. This last property, in particular, 
makes the chameleon field able to satisfy the constraints 
from laboratory tests of the principle of equivalence, be- 
cause in high densities environment, such as terrestrial 
laboratories, the field can be heavy enough so to evade 
them. 

It must be stressed, however, that even with such a 
mechanism, it is very difficult to build a theory with 
a late time cosmology obscrvationally indistinguishable 
from the standard ACDM model. 
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Another interesting consequence of this model conies 
out when considering hnear perturbations of matter and 
their related equation. In the most general case of a 
scalar-tensor theory Q we have: 



5'! + aHSl 



2 



2/32 



1 + 



(8) 



where a is the scale factor, H is the Hubble function, 
6c = S{pme^'^^^'^^') is the matter density contrast and, if 
the field is at the minimum, its mass is = ^00- The 
quantity in brackets can be interpreted as the expression 
for the effective gravitational constant within the context 
of massive scalar field models 1421: 



profile is suppressed by a factor of A. Thus the term 
"thin-shell" comes from the fact that only a portion of 
such a "thin-shell" contributes to the external Yukawa 
profile. A discussion of this issue in the framework of 
/(i?)-gravity is in [l^. 

In summary the main ingredients of our model are: 

• a massive scalar field coupled with ordinary observ- 
able baryonic mass; 

• its mass m, or its interaction length L oc m~^\ 

• its coupling constant with baryonic mass j3 



Geff{a;p,m;k) = 0^1 + 2^^ 



1 



(9) 



In particular the term proportional to results from 
the scalar field-mediated force, which is negligible if the 
physical length scale of the perturbation is much larger 
than the range of the chameleon-mediated force, namely, 
if a/A: ^ m~^. In this case the left hand side of Eq. ([8]) 
is well approximated by 3a^H^Sc/'2 and the matter fluc- 
tuations grow as in general relativity. 

Anyway chameleon theories (i.e. based on Eq. [5]) do 
not behave like linear theories (as anticipated in the in- 
troductory section) of massive scalar fields when massive 
bodies are involved. Varying the action Eq. ([T]) with re- 
spect to the chameleon (scalar) field in a spherically 
symmetric spacetime gives: 



d^ 



2d0 
r dr 



dVr- 



(10) 



because of the effective potential for cj) changes in dif- 
ferent density environments, this differential equation is 
highly non-linear. These non-linear features have been 
studied in [l^, where it was found that for a spherically 
symmetric object of mass Mc and radius Rc surrounded 
by a gas of asymptotic density poo , the profile of the field 
is governed by the so-called "thin-shell" parameter. 



A = 



(11) 



where cp^,^ 



and (/)^j„ are the minima of the effective 
potential outside and inside the object respectively and 
4>c = Mc/SnMpiRc is the Newtonian potential at the 
surface of the object. Thus A is the ratio of the difference 
in (f) potential to the Newtonian potential and quantifies 
how perturbing the object is for the (j> field. If A is large, 
which happens for small objects, then the external (to 
the object) profile of (p is the usual Yukawa profile. 



0(r) 



f3 



AnMpi 



Ale cxp 



(12) 



where niao is the mass of the scalar field. For large 
and compact objects, A is small and the Yukawa 



A. Modified gravitational potential 

Taking the inverse Fourier transform of Eq. Q it is 
straightforward to obtain the corresponding expression 
of the gravitational potential for a point mass distribu- 
tion, ip{r). Remembering that a potential cx in real 
space yields a fc^^ term in Fourier space, we can recog- 
nize in Eq. ^ the point-like gravitational potential per 
unit mass: 



tPir) =-^(l + 2/32e"""-) 



G 



= 1 + 2(3'e 



(13) 



where m is the mass of the scalar field, L cx is the 
interaction range of the modified gravitational potential, 
i.e. the length where the scalar field is effective, and /3 
still being the coupling constant between matter and the 
scalar field. The gravitational potential given in Eq. (fT3|) 
is a point-like one, so that we have to generalize it to ex- 
tended systems if we want to use it for clusters of galaxies 
and LSB galaxies. As we will discuss later, we are going 
to model galaxy clusters as spherically symmetric sys- 
tems: we simply consider the system composed by many 
infinitesimal mass elements each one contributing with a 
point-like gravitational potential. Then, summing up all 
terms, namely integrating them on a spherical volume, 
we obtain a suitable potential. Specifically, we have to 
solve the integral: 



«'(r) 



r'^dr' I sin 9' d9' 



2tt 



duj'i^ir') . (14) 



We make explicit that with an abuse of notation we are 
writing inside the point-like potential r', while it should 
be replaced by |x — x'| = (7-^ + r'^ — 2rr' cos6'')^/^. 

The point-like potential can be split in two terms. The 
Newtonian component is: 



CM 



(15) 



and its extended integral is the well-known expression: 

GM{< r) 



*A,(r) 



(16) 
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where M(< r) is the mass enclosed in a sphere with ra- 
dius r. The correction term coming from the scalar field 



from the integration of the angular part, we have: 

■^c{r)^^2TTG{2(i^L) dr'r'p(r') • 
Jo 



(17) 



e i — e I- 



(18) 



The radial integral is numerically estimated once the 
mass density is given. A fundamental difference between 
such a term and the Newtonian one is that while in the 
latter the matter outside the spherical shell of radius r 
does not contribute to the potential, in the former exter- 
nal matter takes part to the integration procedure, even 
if its contribution is really negligible in most cases. 
At the end, the total potential of the spherical mass dis- 
tribution will be: 



5f(r) = *Ar(r) -f*c(?-) • 



(19) 



As we will show below, for our purpose we need the grav- 
itational potential derivative with respect to the variable 
r; this may not be evaluated analytically so we estimate 
it numerically, once we have given an expression for the 
mass density p(r). While the Newtonian term gives the 
simple expression: 



d$ 



GM{< r) 



dr 



(20) 



the derivative of the corrective potential term is more 
involved. We do not give it explicitly for sake of brevity, 
but only remind that it is an integral-function of the form 



T{r,r') = / dr' f{r,r') 

J Q(r) 



(r) 

from it one has: 

d^(r,r') _ /'^<^^^^,cl/(r,r') 



(21) 



dr 



a(r) 



do: 

- f{r,a{r)) — {r) + 
dr dr 



+ /(r,/3(r))^(r) 



(22) 



Such an expression is numerically derived once the inte- 
gration extremes are given. 

For spiral galaxies, we have the same theoretical appa- 
ratus but a different geometric configuration. Matter in 
spiral galaxies is generally modeled as distributed in a 
thin axis-symmetric disk, so that the extended gravita- 
tional potential is given by: 

/>00 /'27T 

^{r,z)^ R'AR' duj' i;{R',z'). (23) 
Jo Jo 



Even in this case, for being more exact, the couple of vari- 
ables (i?', z') inside the point- like potential should be re- 
placed by |x-x'| = {R'^ + R''^-2RR' cosuj' + {z^z')'^y/^. 
Once the gravitational potential is given, the rotation 
curve for the disk can be easily computed starting from 
the relation M'- 



dv.(i?,.) 



di? 



(24) 



It is important to underline that both in the clusters 
of galaxies case and in the LSB galaxies one, we have to 
perform derivatives with respect of the distance from the 
center of any system, r, of the numerically derived grav- 
itational potential. To be completely rigorous we should 
add a term to all the previously written relations, com- 
ing from the derivative of the function P{p) ^ /3(r). In 
our scalar-field approach we have treated the coupling pa- 
rameter as a constant at all, while in the definition of the 
chameleon mechanism it is also possible and not trivial 
that it may depend on the local density of the gravita- 
tional system one is going to consider. At the same time, 
we do not know what is the possible analytical behavior 
of this quantity or, better: it is one of our purposes trying 
to reconstruct it. It is also evident that if we want to use 
previous relations in the form we have shown before, we 
are implicitly assuming that (3 can satisfy two different 
scenarios: a. it is really constant, and in this case one 
expects not to detect any change in it when comparing 
different gravitational scales; or (as we will shown to be 
our case) b. it is a function of the gravitational scale, but 
its derivative is supposed to be negligible, i.e. d/3/dr sa 0. 



B. Modified distance modulus 

In 0] the Friedmann equation is derived from the ac- 
tion governing the dynamics of the chameleon field (j) in 
the Jordan frame: 



WM., + 1 ^2 ^ y(^) ^ ^ (25) 



with contribution from matter, radiation and the scalar 
field. Making explicit the expression for the Planck mass, 
Eq. ([^5)) becomes: 



= SttG, 



(26) 

We can show that this equation in the chameleon case 
easily converts in the most general expression for a given 
scalar field </>. 

If we assume the chameleon field (j) is in the minimum 
of the effective potential from the earl y s tages of the 
universe on, then we have ^ 1 [42[ during the 

subsequent evolution until today. This also means that 
^p<j>/MpL — \ very high accuracy so it will disap- 
pear from equations and does not have to be considered 
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here. Then, considering that the function F appearing 
in Eq. ([3]) is equal to unity in the scalar field case, we 
also have Gn = G* . At the end the Friedmann equation 
is completely equal to the usual expression: 



3i/^ = SttG 



N 



1 



(27) 



and finally, avoiding radiation contribution while suffix 
Sc corresponds to the scalar field: 



(28) 



where: 



h{z) 



H{z) 



^77), .0 ~ 



Ho 



(29) 



BttG 



N 



Sc,0 



ml 



r + v<j> 



SttGn 

PSc,0 



The function e(z) is unknown; but one knows (assumes) 
that the scalar field works like a cosmological constant on 
cosmological scales, so we may choose it to be constant in 
redshift or one can use more general and extended models 
as the Chevallier-Polarski-Linder (CPL) paramctrization 
[U [5l| usually used to phenomenologically describe dark 
energy fluids. 

Even if we do not have any possibility to discriminate 
between ACDM and a scalar field scenario only by using 
h{z), we have a discriminating tool in the distance mod- 
ulus, the main observable quantity derivable from SNela 
observations, modified from the usual expression by as- 
suming that the gravitational constant can vary with 
time: 



^(z; /3, m; fc) = 5 log [ {1 + z) 
15 



Hz) 
Geff{z;l3,m; k) 



fJ-o 



4 ^°^Ge//(0;/3,m;fc) 



(30) 



In this expression there is an additional term made of 
with the ratio between the value of effective gravitational 
constant at any redshift and the same quantity evalu- 
ated at the present {z = 0). As accurately described in 
(to! , a time- varying gravitational constant can affect light 
curves from SNela by changing both the thermonuclear 
energy release, since the luminosity at the maximum in 
the light curve is proportional to the mass of synthesized 
nickel, and the time scale of stellar explosion. This means 
that by using Eq. ((30)) we are going to test the scalar field 
mechanism on cosmological scales, in particular consid- 
ering the consequent role of a possible variation of the 
effective gravitational constant in the Universe accelera- 
tion rate history. 

In the scalar field theory one has an analytical expres- 
sion for the effective gravitational constant, i.e. Eq. 



which we have modified in the following one for uniform- 
ing all the results: 

Ge//(^;/3,i;A) = G^. (^1 + 2/3^ ^^ii^^^J^) , 

(31) 

which depends on the following variables: 

• the redshift z (a = 1/(1 + z)); 

• the wavelength k (or the length A cx k~^), which 
one could fix or vary on a grid; 

• the intrinsic model parameters, i.e., the coupling 
constant {3 and the interaction length L cx ■m~^, 
which can be constrained with a fitting procedure. 

III. OBSERVATIONAL DATA 

We are going to test the scalar field mechanism on 
three different scale ranges: 

• on cosmological scales, by means of supernovae lu- 
minosity distance; 

• on a Mpc-astrophysical scale, using mass profiles of 
clusters of galaxies; 

• on a kpc-astrophysical scale, analyzing rotation 
curves from spiral galaxies. 

For any of them, we have found out the necessary sample 
data in literature. 



A. Cosmological scale: Supernovae 

SNela are useful because of the possibility to easily 
modify the expression of distance modulus for more gen- 
eral theories with a varying gravitational constant (as in 
the scalar field case). Moreover, the distance modulus is 
the main observable quantity derivable from this kind of 
astrophysical objects. Adding to that, the possibility of 
using data ranging up to redshift values much larger than 
those ones from galaxies or clusters of galaxies (z « 2) 
makes us possible to test and verify a possible temporal 
variation of the gravitational constant (if there is any) 
and so a possible alternative gravity scenario. It is inter- 
esting to underline that with the modified expression of 
the distance modulus we can also verify the coexistence, 
at the same time, of both dark energy and dark mat- 
ter, both explained as different consequences on differ- 
ent scales of the same unified context, namely, the scalar 
field. In fact, in the expression for the distance modu- 
lus we are going to describe in next sections, we employ 
both a term with a dark energy-modeled fluid, coming 
out nothing else that from the effective behavior of the 
scalar field on cosmological scales (which does not dif- 
fer much from the cosmological constant behavior) and a 
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term acting as a dark matter-modeled component, coming 
out from the scalar field working on gravitational scales 
smaller than the cosmological one. In this sense, we have 
also explored how the scalar field process can mimic dark 
matter profiles in clusters and spiral galaxies. 

We use the Constitution sample described in (45l |. 
which is a data set obtained by combining the Union 
data set by [i^ with new 90 nearby objects from the 
CfA3 release described in Q. 

The Union SNela compilation is a data set of low- 
redshift nearby-Hubblc-flow SNela and is built with 
new analysis procedures for working with several het- 
erogeneous SNela compilations. It includes 13 indepen- 
dent sets with SNe from the SCP, High-z Supernovae 
Search (HZSNS) team, Supernovae Legacy Survey and 
ESSENCE Survey, the older data sets, as well as the re- 
cently extended data set of distant supernovae observed 
with HST. After various selection cuts were applied in 
order to create a homogeneous and high-signal-to-noise 
data set, we have final 307 SNela events distributed over 
the redshift interval 0.15 < z < 1.55. 

The CfA3 data set is originally made of 185 multi- 
band optical SNela light curves taken at the F.L. Whip- 
ple Observatory of the Harvard-Smithsonian Center for 
Astrophysics (CfA); 90 of the original 185 objects pass 
the quality cuts of ^3 ^'^id are added to the Union data 
set to form the Constitution one. 

The statistical analysis of Constitution SNela sample 
rests on the definition of the distance modulus given in 
Eq. (|30p . The best fits were obtained by minimizing the 
quantity 

2 / \ raw {fijzj; ^ 0, X,{9i)} - ^iobs{zj))^ 
XsnIMo, a, {till) = 2^ 



C3 



(32) 

where Af = 397 is the number of observed SNela, /i is 
the distance modulus (the observed, fiobs, a-nd the the- 
oretical one, ^{zj; ^0, X, {9i)), cr^ arc the measurement 
variances and {6i} = {/3,L} is the parameters theory 
vector. The nuisance parameter fiQ encodes the Hubble 
parameter and the absolute magnitude M, and has to 
be marginalized over. Giving the heterogeneous origin of 
the Constitution data set, and the procedures described 
in [i^ and [1^ for reducing data, we have worked with an 
alternative version Eq. p2p , which consists in minimizing 
the quantity 



with respect to the other parameters. Here 

AT 

Cl = 



iKzj^f'O = 0, {0i)} - tJ'obs[Zj)f 



E 



(m(^j ; A^O = 0: {^»)} - t^obsiZj)) 



(33) 



(34) 



(35) 



E 



(36) 



It is trivial to see that x%n J^^^ ^ version of XsAfj 
minimized with respect to /io. To that end it suffices to 
notice that 



Xsn(/^o, a, {Oi}) = Cl - 2c2Ato + csAio 



(37) 



which clearly becomes minimum for /^o = C2/C3, and so 
we can see = xiN(/^o = 0, A, {Oi\). Furthermore, one 
can check that the difference between x|jv ''^^'^ x|jv 
negligible. 

We minimize the x-square using the Markov Chains 
Monte Carlo Method (MCMC) and testing their conver- 
gence with the method described by [13]. The ia confi- 
dence levels are easily estimated deriving them from the 
final samples, using the 15.87-th and 84.13-th quartiles 
(which define the 68% confidence interval) for i = 1, the 
2.28-th and 97. 72-th quartiles (which define the 95% con- 
fidence interval) for i = 2 and the 0.13-th and 99.87-th 
quantiles (which define the 99% confidence interval) for 
i = 3. 



B. Galaxy Cluster Sample 

Clusters of galaxies are uniquely useful tracer of cos- 
mological evolution and so ineludible tests in the field 
of alternative gravities other than general relativity [9l| . 
They are fundamental tracers for two main features. 
First of all, they are the largest gravitational objects 
whose masses can be adequately measured, and the 
largest objects to have undergone gravitational relax- 
ation and entered into virial equilibrium. Second, clus- 
ters are essentially "closed boxes" that retain all their 
gaseous mass content because their gravitational wells 
are much deep. The most accepted paradigm is that 
clusters of galaxies are mostly made of coUisionless cold 
dark matter particles (CDM model) and arc virializcd 
systems from scale-free Gaussian initial density pertur- 
bations. The CDM paradigm and the numerical simula- 
tions make clear predictions for the structure of clusters 
of galaxies; comparisons of these predictions with the re- 
sults of high-quality observations is a necessary consistent 
check and any significant deviation can place important 
constraints on their theoretical model but also on cosmo- 
logical models and, as in our case, on the exploration of 
different gravity theories. 

The formalism described in § III Al can be applied to 
a sample of 12 galaxy clusters. We shall use the clus- 
ter sample studied in [s^ which consists of 13 low- 
redshift clusters spanning a temperature range 0.7 -^- 
9.0 keV derived from high quality Chandra archival data. 
In all these clusters, the surface brightness and the gas 
temperature profiles are measured out to large radii, so 
that mass estimates can be extended up to r5oo or be- 
yond. 
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Clusters of galaxies are generally considered self-bound 
gravitational systems with spherical symmetry and in hy- 
drostatic equilibrium if virialized. The last two hypothe- 
ses are still widely used, despite of the fact that it has 
been widely proved that most clusters show more com- 
plex morphologies and/or signs of strong interactions or 
dynamical activity, especially in their innermost regions 

Under the hypothesis of spherical symmetry in hydro- 
static equilibrium, the structure equation can be derived 
from the collisionless Boltzmann equation 



2/3gas(^) / 2 
Z (.^r- 



-Pgasir) 



(38) 



where 4" is the gravitational potential of the cluster, ar 
and ag^uj are the mass- weighted velocity dispersions in 
the radial and tangential directions respectively, and p is 
the gas-mass density. For an isotropic system, it is 

o-r = cr0,oj ; (39) 
the pressure profile can be related to these quantities by 



P(r) 



Substituting Eqs. 1^ - (gO]) into Eq. 
an isotropic sphere. 



dP(r) 



-Pgasir) 



d*(r) 



(40) 

we have, for 
(41) 



For a gas sphere with temperature profile T(r), the ve- 
locity dispersion becomes 



fcT(r) 

fJ-Ulp 



(42) 



where k is the Boltzmann constant, fj, w 0.609 is the mean 
mass particle and nip is the proton mass. Substituting 
Eqs. dini) - dH]) into Eq. (gl]), we obtain 

d fkTir) \ d^ir) 

Pgasir) = -pgasir)- 



dr \ /imp 
or, equivalently, 

dl'(r) _ fcr(r 



dr 



firripT 



dr 



din Pgasir) dlnr(r) 



(43) 



d In r d In r 

Now the total gravitational potential of the cluster is: 

*(r) = *Ar(r)+*c(0- (44) 

It is worth underlining that if we consider only the stan- 
dard Newtonian potential and its derivative in Eq. (|43p. 
the total cluster mass Md^Nir) (the standard estimation 
of clusters mass in a CDM scenario) is composed by gas 



mass -|- galaxies mass -I- dark matter and it is given by 
the expression: 



Ma^Nir) = Mgasir) + Mgaiir) + Mouir) 



kTir) 
n'' 



din Pgasir) d In T(r) 



dlnr 



d In r 



(45) 



Generally the galaxy contribution is considered negligi- 
ble with respect to the other two components so we have: 

Md,jv(r) « Mgasir) + MoKiir) « 



kTir) 
pirip 



din Pgasir) dlnr(r) 



dlnr 



d In r 



Since the gas-mass estimates are provided by X-ray ob- 
servations, the equilibrium equation can be used to derive 
the amount of dark matter and to reconstruct its spatial 
distribution in a cluster of galaxies. 

Inserting the previously defined extended- corrected po- 
tential of Eq. dMl into Eq. (gS]), we obtain: 



d^- 



N 



dr 



d*c 
dr 



kTir) 
fj.mpr 



din Pgasir) dlnr(r) 



din) 



din? 



(46) 

from which the extended- corrected mass estimate fol- 
lows: 

r^ d'i'cir) 



cLEC 



ir) 



G dr 
kTir) 
prUpG 



(47) 



din Pgasir) dlnT(r 



din) 



dhu 



Since the use of a corrected potential avoids, in princi- 
ple, the additional requirement of dark matter, the total 
cluster mass, in this case, is given by only the baryonic 
matter counterparts: 



Md^Ecir) = Mgasir) + Mgaiir) , 



(48) 



that can be entirely evaluated by observational data. The 
mass density in the ^'c term is 



Pcl,Ecir) = Pgasir) + Pgalir) , 



(49) 



with the density components derived from observations. 
Considering that the right term in Eq. (|T7|) is the total 
Newtonian mass estimation for a cluster of galaxies, we 
easily derive that the corrective term in the gravitational 
potential works in mimicking an effective dark matter 
contribution: 



r^ d^fc 
G~dr 



ir) = Md^Nir) - Md^Ecir) 



(50) 



But in our approach, instead of requiring new kind of 
particles, it arises by the interaction of baryonic matter 
with the scalar field scalar field. 

We have hence performed a best-fit analysis of the the- 
oretical estimation of dark matter, Eq. ([50)1 



Mdm,thi-r;P,L) = Meffir;P,L) 



r^ d4'c(r) 
G dr 



ir) (51) 
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which depends on scalar field parameters through the po- 
tential ^c", versus the same but observationally-derived 
quantity, 



Mdm,obsir) = McLNir) - McLEcir) 



(52) 



We underline here that we could not fit directly the ob- 
served baryonic mass because of the great difference in or- 
der of magnitude between Md^Nif) a-nd Mci^ec{''')i work- 
ing this last one like a small perturbation to the total mass 
estimation. It is clear that the term corresponding to our 
theoretical derived dark matter quantity is much bigger 
than the baryonic contribution, and even a small and 
acceptable deviation of only 1% in it would have been 
translated in a larger deviation of 10% in the baryonic 
one. 

Since not all the data involved in the above estimations 
have measurable errors, wc cannot perform an exact x~ 
square minimization. Actually, wc can minimize the 
quantity: 



1 



(53) 

where M is the number of data and Up — 2 is the free 
parameters number of the model and {6i\ = As 
usual wc find the minimum in x-square running MCMCs. 
Even if the convergence is achieved after few thousand 
steps of the chain, we have decided to run longer chains 
of 10^ points to reduce the noise from the histograms 
and avoid under- or over- estimations of errors on the 
parameters. 



1. Gas Density Model 

The gas density distribution of the clusters in the sam- 
ple is described by the analytic model proposed in [90l |. 
Such a model modifies the classical /3— model to represent 
the characteristic properties of the observed X-ray sur- 
face brightness profiles, i.e. the power-law-type cusps of 
gas density in the cluster center, instead of a flat core and 
the steepening of the brightness profiles at large radii. 
Eventually, a second /3— model, with a small core radius, 
is added to improve the model close to the cluster cores. 
The analytical form for the particle emission is given by: 



(r/re)-" 



1 



(l + r2/r2)3/3-"/2 [l + r-y/riyh 



'02 



(l + rVr22)3/32 ' 



(54) 



which can be easily converted to a mass density using 
the relation: 



1.4 

Pgas = riT ■ fJ-rrip = — rieTOp 



(55) 



where nx is the total number density of particles in the 
gas. The resulting model has a large number of parame- 



ters, some of which do not have a direct physical interpre- 
tation. While this can often be inappropriate and com- 
putationally inconvenient, it suits well our case, where 
the main requirement is a detailed qualitative descrip- 
tion^of the cluster profiles. 

In ^9Qj, Eq. ([M)) is applied to a restricted range of dis- 
tance from the cluster center, i.e. between an inner cut- 
off rmim chosen to exclude the central temperature bin 
(w lO-f-20 kpc) where the ICM is likely to be multi-phase, 
and rdet, where the X-ray surface brightness is at least 
3(7 significant. We have extrapolated the above function 
to values outside this restricted range using the following 
criteria: 

• for r < r„ii„, we have performed a linear extrapo- 
lation of the first three terms out to r = kpc; 

• for r > rdet, we have performed a linear extrapo- 
lation of the last three terms out to a distance f 
for which Pgasij) = Pc, Pc being the critical den- 
sity of the Universe at the cluster redshift: pc ~ 
Pc,o ■ (1 + z)^- For radii larger than f, the gas den- 
sity is assumed constant at pgas{r). 

We point out that, in Table |lTl the radius limit is al- 
most the same as given in the previous definition. When 
the value given by [9^ is less than the cD-galaxy radius, 
which is defined in the next section, we choose this last 
one as the lower limit. On the contrary, Tmax is quite 
different from rdet- it is fixed by considering the higher 
value of temperature profile and not by imaging meth- 
ods. 

We then compute the gas mass AIgas{r) and the total 
mass Mci^Nir), respectively, for all clusters in our sam- 
ple, substituting Eq. into Eqs. and re- 
spectively; the gas temperature profile is described in 
details in § IIIIB2I The resulting mass values, estimated 
at r = rmax , are listed in Table HT] 



2. Temperature Profiles 



As stressed in § IIII B 11 for the purpose of this work, 
we need an accurate qualitative description of the radial 
behavior of the gas properties. Standard isothermal or 
polytropic models, or even the more complex one pro- 
posed in [g^l, do not provide a good description of the 
data at all radii and for all clusters in the present sample. 
We hence describe the gas temperature profiles using the 
straightforward X-ray spectral analysis results, without 
the introduction of any analytic model. 
X-ray spectral values have been provided by A. Vikhlinin 
(private communication). A detailed description of the 
relative spectral analysis can be found in [89j . 



3. Galaxy Distribution Model 

The galaxy density can be modelled as proposed by 
[3| . Even if the galaxy distribution is a pomt-distribution 
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instead of a continuous function, assuming that galaxies 
are in equilibrium with gas, we can use a /3— model, oc 
r~^, for r < Rc from the cluster center, and a steeper 
one, cx r~^-^, for r > Rc, where Rc is the cluster core 
radius (its value is taken from [1^). Its final expression 
is: 



Pgaiir) 









Pgal,l ■ 


1 + 


(*)1 








Pgaia ■ 


1 + 


(*)^: 



r<Rc 

(56) 

r>Rc 



where the constants Pgai,i and Pgai,2 are chosen in the 
following way: 

• provides the central number density of galaxies 
in rich compact clusters for galaxies located within 
a 1.5 h~^Mpc radius from the cluster center and 
brighter than 771,3 + 2™ (where 7713 is the magni- 
tude of the third brightest galaxy): Ugai^o ~ IQ^h^ 
galaxies Mpc~'^. Then we fix Pgai,i in the range 
- 1034 ^ 1036 kg/kpc3. 

• the constant PgaL2 has been fixed with the only 
requirement that the galaxy density function has 
to be continuous at Rc- 

For any cluster we assume that the galaxy population 
also consists in a cD galaxy, a giant elliptical galaxy with 
a diffuse envelope which is generally located at the center 
of clusters and whose typical mass is in the range 10^^ 
IQ^^Mq. The cD galaxy density has been modeled as 
described in [t^ ; they use a Jaffe model of the form: 



PCDgal 



Po,J 



1 



(57) 



where rc is the core radius while the central density is ob- 
4 

tained from Mj = -ttRIpq^.j. The mass of the cD galaxy 

O 

has been fixed at 1.14x10^2 Mq, with = i?e/0.76, with 
Rc = 25 kpc being the effective radius of the galaxy. The 
central galaxy for each cluster in the sample is assumed 
to have approximately this stellar mass. 

We have tested the effect of varying galaxy density with 
the central density parameter Pgai,i in the above range 
~ lO^"* -;- 10^6 kg/kpc3 on the cluster with the lowest 
mass, namely A262. In this case, we would expect greater 
variations with respect to other clusters; the result is that 
the contribution due to galaxies and cD-galaxy gives a 
variation < 1% to the final estimate of fit parameters. 

Finally, we have assumed that the total galaxy- 
component mass (galaxies plus cD-galaxy masses) is w 
20 25% of the gas mass: in [t^I, the mean fraction of 
gas versus the total mass (with dark matter) for a cluster 
is estimated to be 15 20%, while the same quantity for 
galaxies is 3 5%. This means that the relative mean 
mass ratio gal-to-gas in a cluster is ss 20-;- 25%. We have 
varied the parameters Pgai,i, Pgai,2 and Mj in their pre- 
vious defined ranges to obtain a mass ratio between total 



galaxy mass and total gas mass which lies in this range. 
At the end the cD-galaxy is dominant with respect to the 
other galaxies only in the inner region (below 100 kpc). 
As already stated in § IIIIB 11 cluster innermost regions 
have been excluded from our analysis and so the contri- 
bution due to the cD-galaxy is practically negligible in 
our analysis. The gas is, as a consequence, the dominant 
visible component, starting from innermost regions out 
to large radii, being galaxy mass only 20 25% of gas 
mass. 



4- Uncertainties on mass profiles 

Uncertainties on the cluster total mass profiles have 
been estimated performing Monte-Carlo simulations (63j . 
We proceed to simulate temperature profiles and choose 
random radius-temperature values couples for each bin 
which we have in our temperature data given by j89| . 
Random temperature values have been extracted from 
a Gaussian distribution centered on the spectral values, 
and with a dispersion fixed to its 68% confidence level. 
For the radius, we choose a random value inside each 
bin. We have performed 2000 simulations for each clus- 
ter and perform two cuts on the simulated profile. First, 
we exclude those profiles that give an unphysical negative 
estimate of the mass: this is possible when our simulated 
couples of quantities give rise to too high temperature- 
gradient. After this cut, we have « 1500 simulations for 
any cluster. Then we have ordered the resulting mass 
values for increasing radius values. Extreme mass esti- 
mates (outside the 10-^90% range) are excluded from the 
obtained distribution, in order to avoid other high mass 
gradients which give rise to masses too different from real 
data. The resulting limits provide the errors on the to- 
tal mass. Uncertainties on the electron-density profiles 
have not been included in the simulations, being them 
negligible with respect to those of the gas-temperature 
profiles. 



C. Low surface brightness galaxies 

For the analysis of galactic scales we have used a sam- 
ple of the so-called low surface brightness (LSB) galaxies 
and dwarf galaxies. It is still unclear to what extent 
rotation curves in bright spiral galaxies may give clues 
about the profile of both visual and dark matter, mainly 
because they are poor in gas content, so that rotation 
curves can be hardly detected out to sufficiently large 
radii where they are supposed to be dark matter dom- 
inated. Moreover they also show some typical complex 
features, such as extended spiral arms or barred struc- 
ture that can lead to consistent non-circular motions and 
thus making difficult the interpretation of data. On the 
other side, LSB galaxies are supposed to be dark mat- 
ter dominated at all radii, and therefore the analysis of 
their rotation curves can yield important clues about it. 
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Effectively, LSB galaxies exhibit a large discrepancy be- 
tween the detectable and the Newtonian dynamics mass 
within the optical disk [H [13, [13, HI . They are also a 
little challenge in the framework of CDM model because 
predictions from CDM based simulations have revealed 
disagreement with observational profiles of several dwarf 
galaxies [55|. In particular, data indicate much less cuspy 
distributions of matter than the simulations, and possi- 
ble solutions such as feedback effect due to star formation 
have been excluded [s^ by the observed low star forma- 
tion ratio, at present and in the past, from which it can be 
deduced that star formation rate has never been impor- 
tant enough in these galaxies to modify their structure. 
Then, alternative solutions to Newtonian dynamics are 
also possible. 

The chosen galaxies for our analysis come from a sam- 
ple of 15 elements with high resolution Ha/HI rotation 
curves extracted from the lar ger sample described in jsij . 
This sample was selected by [l^l using the criteria of the 
contemporary availability of data on rotation curves, on 
surface photometry in the R band and on surface gas- 
mass density. Rotation curves were derived from spec- 
trographic observations performed by the authors them- 
selves of [m , while the photometry of the stellar disk and 
Ha/HI surface densities are collected from literature. For 
a more complete and detailed discussion one can mainly 
refers to 

The final sample does not constitute a complete sam- 
ple of dwarf and LSB galaxies, but contains galaxies in a 
wide range in luminosity and surface brightness for which 
high-quality rotation curves are available. Therefore it is 
well suited for testing scalar field mechanism in late- type 
dwarf and LSB galaxies. 

Moving to the way of modeling our spiral galaxies, we 
have to specify the properties of stellar and gas distribu- 
tion. Stars are generally supposed being distributed in a 
thin and circularly symmetric disk, with the surface den- 
sity derived from the observed surface brightness 
distribution through the relation: 



S](i?) = 10 



-O.4(M(fl)-AiH,0-C) 



(58) 



where (ir^q = 4.46 is the solar magnitude in the R band 
and C = 21.572 is the needed constant for converting 
surface brightness from magnitude units, mag/arcsec^, 
to linear units, Lq/ pc^. The luminosity surface de nsity 
is commonly fitted with the exponential disk model [4l| : 



S](i?) = I]oCxp(-i?/i?rf) , 



(59) 



where Sq is the central surface brightness and Rd is the 
disk scale length. Of course, in determining our theo- 
retical rotation curve, we need the stellar mass density, 
that can be obtained from the luminosity one by simply 
multiplying it for the stellar mass-to-luminosity ratio, Y*, 
which is, together with our scalar field parameters, the 
third free parameter having to be constrained with the 
fitting procedure. 

Modeling the gas density is more complicated because 
we do not have an analytical function able to describe 



its behavior at all radii and because the profile is very 
disturbed. Using [7§| plots and images, we fit the outer 
radii profile with a linear relation (in magnitude unities), 
while the inner one is generally fitted by simply inter- 
polating data, with any analytical expression (generally 
a polynomial one) which well fits data points. We then 
check if the model works well by comparing the obtained 
total gas mass with the same quantity but evaluated by 
extrapolation from observational data. We verified that 
we only need very small normalization constants, in the 
range (0.95, 1.05), to fully match data values. Finally, we 
multiply the gas density with the factor 1.4 for including 
helium contribution. 

With these model components, the general expression 
for the rotation velocity 



vliR)^R 



d*(i?,z) 



di? 



2 = 



can be decomposed in the following terms: 
vliR) = vlj,{R)+vlciR) , 



with 



and 



(R) = R 



d*^(i?,z) 



(i?) = R 



di? 



dR 



(60) 



(61) 



(62) 



(63) 



the Newtonian and the corrective contributions to the 
rotational velocity from the respective terms in the total 
gravitational potential. Anyone of the previous terms can 
then be decomposed in two different component elements, 
concerning the two mass components, namely stars and 
gas. For stars we have: 



star 2 



^d*?*"''(i?,z;n 
il — 



di? 



(64) 



2=0 



where we make explicit that this term depends on the 
free parameter which appears in the stellar density 
inside the potential expression. When this term is not 
present in the available data, we use the expression given 
in for the case of an exponential disk and with 1", = 1 , 
i.e. 



star 2 



(i?) =4^GSoi?d(y)2 



(65) 



with y = R/2Rd. For gas we have 



gas 2 



(i?) = i? ■ 



dvf-rc^,^) 



di? 



(66) 



and, as in the stellar case, when this kind of data is not 
available, we derive it in a numerical way by using the 
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modeled gas density. From the corrective term to the 
potential we have 



d*^*'^'-(i?,z;K,/3,X) 



dR 



2=0 



and 



qas 2 



{R-P,L)^R 



dv|/^''^(i?,z;/3,i) 



di? 



(67) 



(68) 



both these two quantities are derived numerically. Fi- 
nally, the total rotation velocity is the sum in quadrature 
of all these elements, i.e. 



and the chi-square function is; 

A/" 



xLB(m)-E 



(wM/.(-R„{^«)}-«c,ofcs(i?0)' 



(69) 



(70) 



where M is the number of data, af are the measurement 
variances and {0^} — {/3,Z/,Y^,} is the parameters the- 
ory vector. Even in this case, like for supernovae and 
clusters, we use MCMC method for minimizing the chi- 
square function and deriving errors on the fit parameters. 



IV. RESULTS AND DISCUSSION 

In this section we are going to discuss the obtained 
results with the goal to achieve a comprehensive back- 
ground where scalar field works at various scales. In par- 
ticular, we are going to search for possible trends and 
correlations between observable quantities and theoreti- 
cal scalar field parameters in the various class of the con- 
sidered astrophysical objects. SNela will be discussed 
alone because of the difficulties in interpreting their re- 
sults in this general context. On the other hand, clusters 
and galaxies seem to be tied in a common picture by a 
rescaling process. 



A. Supernovae: results 

Difficulties come from the SNela analysis. In this case 
we have performed MCMC analysis leaving the scalar 
field parameters free, with the only minimal requirement 
that L > 0, being it a length, and that /? > 0, because of 
only the term 0^ appears to be involved in the Eq. ([51]) 
and no possibility to discriminate between positive and 
negative values is given. 

The parameter A may be considered as a length related 
to stellar scale, to typical supernovae explosion, and/or 
to stellar formation. Since its exact value is unknown, 
we have varied it on a grid ranging from A = 10~^h~"'^ as 



the minimum, and corresponding to a length of ~ 1 kpc 
(assuming Hq = 74.2 as in [HO]), up to A = Ih^^ as the 
maximum, and corresponding to a length « 1 Mpc. 

In Table |T] we have also reproduced the values of chi- 
square function evaluated at the best fit points: they are 
completely equal to the same quantity evaluated for a 
CPL model with rim.Oi ^-nd Wa as fit parameters and a 
gravitational constant fixed at its usual Newtonian value 
[constant in time and scale). For CPL model we have 
obtained the following results: 



nr. 



0.273: 



-0.091 
-0.133 



Wo 



U.yo^_0 275 

- 003+0-201 



Xbest = 465.665 



(71) 



If we want to compare the reliability of our model against 
the CPL one, we can use a Bayesian type test as BIC. It 
is defined as BIC = -21n/: + fclniV with -21n£ be- 
ing the chi-square value and £ the likelihood function, k 
the number of parameters of the model and N the num- 
ber of points in the dataset. The best model has the 
lowest BIC; in particular, when comparing two different 
models, if the difference between the two BIC values is 
ABIC < 2, than there is no significant difference between 
the models; if 2 < ABIC < 5, this difference is substan- 
tial; if 5 < ABIC < 10 there is a "strong" evidence 
in favor of the model with the lowest BIC value; while 
for ABIC > 10 this evidence is "decisive" (following the 
most used value scale, the "Jeffreys' scale", defined and 
reported in [47j). 

If we consider that in the CPL case we have three fit 
parameters as in our scalar field model, we can easily de- 
duce that there are no significative differences between 
the CPL and the scalar field model. 
Things change slightly if we consider a ACDM model 
with a Newtonian gravitational constant and only one 
free parameter, fi,n.o- In this case, even if the chi-square 
has the same value of our scalar field model, we have only 
one parameter, so that ACDM model with a Newtonian 
gravitational constant is strongly favored with respect to 
a scalar field model with a varying gravitational constant 
(we may underline that this is an obvious consequence of 
the BIC definition, where the fclnTV tends to support 
models with a smaller number of parameters). We will 
return later on the chi-square values for a further discus- 
sion. 

Now we can take a look to the best fit values of model 
parameters. For what concerning the matter content pa- 
rameter, we can observe that for fim.o they are slightly 
higher than the latest estimation [5^, Vlm^o ~ 0.24, from 
combining WMAP results with BAO and and the Hub- 
ble constant measurements, but they are aligned with 
usual values derived from only-SN analysis, which gener- 
ally lead to higher values for this parameter [l5| . 
We have also to stress that such a high value for rim.o, 
although comparable with the usual contribution gener- 
ally attributed to dark matter, is not contradictory at 
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Table I. Supernovae. Column 1: A value. Column 2: chi square value evaluated at the best fit values for fitting parameters. 
Column 3: present matter content, ^lm,o- Column 4: coupling parameter /3 fi-om scalar field (la confidence interval). Column 
5: gravitational length L from scalar field (la confidence interval). 

^ Xbeat ^m,0 P L 

[h'^) (Mpc) 



10-3(«lfcpc) 465.610 0.292i;j;°g 0.135t°;*|^ 0.145115;?^^ 

10~^ (f» 10 fcpc) 465.659 0.295t!5:o2l Q.12Qto'f^ 0.1581^;^^ 
10-^ 100 fcpc) 465.633 O.SOllgigje 

10° (^1 Mpc) 465.707 Q.2Q?,t°o°^\ 0.129t^;g9;^ 0.159 



*-0.083 '-'•^'^-0.136 

1+0.991 
-0.123 



all with our intent of explaining dark matter on smaller 
than cosmological scales and dark energy on cosmologi- 
cal scales with the same source, namely, the scalar field 
effective scalar field. 

In fact we can easily consider that there are two elements 
contributing to the rim.o value: one coming from ordinary 
baryonic matter, fi^.o ~ 0.04, and one coming from the 
scalar field, thus acting as a dark matter-type component 
and scaling as (1 + zj^ . As the discussion in next pages 
will show, this is a feasible possibility. Then, we have a 
contribution from scalar field to the acceleration of uni- 
verse through an effective dark energy component, for an 
amount equal to the usually derived one (fic/i.o ~ 0.70) 
and acting as a cosmological constant. 
Discussions about the scalar field parameters arc more 
tricky: they do not show significative statistical changes 
while varying A. By applying Eq. pop to the Constitu- 
tion data set we are implicitly assuming that the A length 
is the same for all supernovae; if we assume that these 
objects constitute an homogeneous astrophysical family, 
our hypothesis looks general and quite good. 
The parameter /? quantifies the coupling of the scalar field 
with ordinary matter and mainly measures how much the 
effective gravitational constant deviates from the usual 
Newtonian one. From the cosmological analysis it is in 
the interval [0.114,0.135]. On the other side, the length 
L may be considered in this case as the "minimal" length 
for which the variation of the gravitational constant has 
"cosmological" effects ( i.e. detectable with the Hubble 
SNela diagram) and with the scalar-field mimicking a 
dark energy-type component. Our analysis shows that 
such a scale is ^ 100 kpc. 

Anyway, we have to combine all these results with the 
evidence that chi-square values are completely equal for 
any A, then we have also to consider the possibility that 
our analysis is not well-based. In the top panel of Fig. ^ 
we have plotted how much the modified distance modulus 
in Eq. pOI) differs from the usual expression. 



^(z) = 51og (1 + 2) 



+ MO • (72) 



/o Kz), 

If we consider that the modified distance modulus can be 



written as: 

fi{z;^,L;\) = fi{z) 



15 . Ge//(z;/3,^;A) 
4 °^Ge//(0;/3,i;A) 



(73) 



the correction coming from the effective scalar field grav- 
itational constant is negligible when compared to the 
usual distance modulus expression, being at most 
0.08%. Facing with this question, we have two possi- 
bilities: 

1. our results do not come from an effective test of 
a scalar field mechanism, but are got stuck by the 
impossibility to detect changes in the gravitational 
constant, being the distance modulus mainly driven 
by the flm.o parameter in Eq. (j72p : 

2. we arc effectively testing a scalar field mechanism, 
but then wc have the problem of including the ob- 
tained values for scalar field length L in a consistent 
theoretical background. 

As we can see in the middle and bottom panels of 
Fig. ([2), we have different scenarios depending on the 
value of A: 

• for A = 10""^ wc have an effective gravitational 
constant which is practically constant, exhibiting 
a change of only ^ 0.003% in the redshift range 
[0.005, 1] and assessing at a value which differs 
from the Newtonian gravitational constant value 
for ^ 3.6%. In this case the relative contribution 
of the corrective term with respect of the usual one 
is « 10-5%; 

• for A = 10^^, we have again a practically constant 
effective gravitational constant, with a change of 
^ 0.7% in the SNela redshift range, and assessing 
at a value which differs from the Newtonian grav- 
itational constant value for ~ 2.85%. In this case 
the relative contribution of the corrective term with 
respect of the usual one is really small, sa 0.0006%; 

• for A = 10^^ we have a more sensible rise of effec- 
tive gravitational constant in the supernovae red- 
shift range, reaching an asymptotic value bigger 
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Figure 2. Continuous line is for A = 10"^ (/3 = 0.135; L = 0.145 Mpc); dashed line is for A = 10"^ (/3 = 0.120; L = 0.158 
Mpc); dot-dashed line is for A = 10"^ (/3 = 0.114; L = 0.174 Mpc); dotted line is for A = 10° (/? = 0.129; L = 0.159 Mpc). 
The vertical dotted line shows the maximum redshift from SNela sample. Top panel. Deviations of the corrected version 
of distance modulus from the usual one. Middle panel. Ratio between the effective gravitational constant, Eq. ((9|, and the 
Newton gravitational constant value. This ratio quantifies the temporal evolution of the deviation of scalar field mechanism 
from Newtonian gravity. Bottom panel. Ratio between the effective gravitational constant, Eq. 0, and the present value of 
the same quantity. This ratio quantifies the temporal evolution of the effective gravitational constant. 
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than Newtonian one for ~ 2.6% at z > 50. Follow- 
ing the large change in the redsfhift dependence of 
the effective gravitational constant, which can vary 
for 1% even in the supernovae redshift interval, the 
ratio between the corrective and the usual expres- 
sion of distance modulus in this case is « 0.028%; 

• for A 1 the rising is less pronounced in the super- 
novae redshift range with respect of previous case, 
but goes over at z 15 and reaches its asymptotic 
value at z ^ 500, being it ^ 3.32% bigger than 
Newtonian gravitational constant. In this case the 
corrective distance modulus term reaches a maxi- 
mum deviation of « 0.08% around z ^ 50; in the 
same range the effective gravitational constant can 
vary for « 3%. 



B. Clusters of galaxies: results 

When considering clusters of galaxies we remind that 
in this case we left free the scalar field parameters, /? and 
i, with only the minimum requirement of their positive- 
ness. 

As it is possible to sec by simple visual inspection, the 
only bad fit corresponds to the cluster RXJ1159: using 
the modelled matter densities described in § piIB \\ - 
pil B 3[) . we obtain a too fast decreasing mass profile in 
the inner region reaching unphysical negative values. For 
this reason we will not consider it anymore in our con- 
siderations. 

On the contrary, for all the other clusters we have good 
results, with mass estimations corresponding in the Ict 
confidence level. Errors contours reported in Fig. ([3]) - 
have two contributions: the main one comes from the 
statistically derived errors on mass observations, as de- 
scribed in § pil B 4p and which produce the larger and 
irregular borders of the Icr confidence level; the smallest 
one comes from errors on fitting parameters. 

What can be first noted is that for lower scales there is 
a large deviation between our theoretical estimation and 
the observed one. Typically one can locate this break- 
point in the range [100, 150] kpc. This is not a new thing 
when describing clusters of galaxies with modified grav- 
ities and it is not an intrinsic failure of our theoretical 
model. Similar issues are present in [2l|, where f{R)- 
gravity models are applied to clusters of galaxies. The 
same situation is in [Ij]: they use the the Metric - Skew 
- Tensor - Gravity (MSTG) as a generalization of the 
Einstein General Relativity and derive the gas mass pro- 
file of a sample of clusters with gas being the only bary- 
onic component of the clusters. They consider some clus- 
ters included in our sample (in particular, A133, A262, 
A478, A1413, A1795, A2029, MKW4) and the same fail- 
ing trend is found for r < 200 kpc: they overestimate gas 
mass in the inner regions with respect of the expected 
estimation from X-ray observations. In the same work 
there is also an interesting note about MOND theory ap- 



plied to clusters of galaxies: even if it is not possible to 
assess it really fails, surely we can see that MOND in 
clusters does not solve the dark matter problem because 
it again requires including a mass contribution other than 
the observed one. 

The reason for this different behavior in the inner re- 
gions is in the break of the hypothesis of hydrostatic equi- 
librium. If the hypothesis of hydrostatic equilibrium is 
not correct, then we are in a regime where the funda- 
mental relations Eqs. ([55)) - are not working. As 
discussed in [sl], the central - 70 kpc - region of most 
clusters is strongly affected by radiative cooling and thus 
its physical properties cannot directly be related to the 
depth of the cluster potential well. This means that, in 
this region, the gas is not in hydrostatic equilibrium but 
in a multi-phase state. In this case, the gas tempera- 
ture cannot be used as a good standard tracer. Among 
the main phenomena which causes this we have cooling 
flows, merger effects and asymmetric shapes. In particu- 
lar, cooling flows produce a decrease in the temperature 
profile and consequently local higher gas densities which 
cannot be related directly to gravitational effects. 
A coherent behavior is shown in our plots. Fig. ([3]) - 
(|31)- We remind that there the distribution of dark mat- 
ter is represented; higher densities from cooling flows 
produce higher non-gravitational-produced gas mass pro- 
flies, which result in a decrease of the dark matter one. 
In our case, we have that our theoretical dark matter 
profile (i.e. the effective dark matter mimicked by the 
different coupling of scalar field with baryonic mass) is 
higher than the observationally derived one (see for ex- 
ample AbeU 262, Abefl 383, AbeU 478, Abefl 1413, AbeU 
1991, AbeU 2029; while for AbeU A133 and AbeU 1795 
one can perceive the same trend but unfortunately the 
data do not extend too small enough radii). 
However, a more detailed modelling for inner regions is 
out the purpose of this work, while we are here inter- 
ested to show that the scalar field mechanism can be a 
valid alternative to dark matter in order to explain clus- 
ter dynamics. In this sense it is very illuminating that 
on the most part of objects in our sample there is a very 
good agreement between the scalar field model and the 
observationally derived dark matter profiles on a wide 
range of distances from the centers for any cluster, ap- 
proximately in the interval [100; 1000] kpc, or at least up 
to the maximum radii coming from observations. 

Giving a more detailed glance to the absolute values of 
scalar field parameters for clusters, we can see that the 
coupling constant /3 looks very well constrained in a re- 
ally narrow range, [1.888; 3.259]; while the gravitational 
length L can vary in the interval [148.931; 1487.800] kpc, 
which seems consistent with the range depicted by other 
used characteristic clusters scale, as the virial radius in a 
ACDM context. 

To be more precise, if we consider (and exclude) three 
peculiar cases, these intervals can even be constrained 
in narrower windows: for /3, the interval may be 
[2.047; 2.786], while for L it may be [147.977; 1323.890]. 
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Figure 3. Dark matter profile vs radii for clusters of galaxies. Dashed line is the observationally derived estimation of dark 
matter, Eq. (|52p : solid line is the theoretical estimation for the effective dark matter component, Eq. 1)51^: dot-dashed lines are 
the 1-a confidence levels given by errors on fitting parameters plus statistical errors on mass profiles as discussed in § IIII B 41 



In [iil, MKW4 and RXJ1159 arc recognized as the 
strongest outliers in temperature profile, exhibiting very 
compact cooling regions and having a temperature peak 
smaller than other sample objects and located in inner 
region at r w 50 kpc. RXJ1159 is better classified as an 
X-ray Over-Luminous Elliptical Galaxy; optically, this 
object appears as a nearly isolated elliptical galaxy but 
its X-ray luminosity and extent is typical of poor clus- 
ters. MKW4 is considered a group of galaxies or a poor 
cluster too, so we have to consider them as different 
object from the rest of the sample. Moreover we cannot 
give any conclusion about this gravitational class (group 
of galaxies) because only one object has no statistical 
weight for any analysis. 



Always in , the cluster Abell 2390 seems to elude typ- 
ical clusters scaling relations; this is due to its unusual 
central cool region, which extends up to r ~ 400 kpc, 
probably because the cold gas is pushed out from the 
center by radio lobes. 

From these considerations, it is interesting to note that 
just these three objects result having the most different 
values for the coupling constant /3 with respect of the 
other clusters and even seem to exhibit a peculiar trend. 
MKW4 (and RXJ1159) have more compact inner cool 
regions and this is associated with: 1. higher values for 
/?, which means a larger coupling of scalar field with or- 
dinary matter (/3 acts like a concentration parameter) 
and 2. smaller value for the gravitational length. On the 
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Figure 4. Same of Fig. (|3]). 



contrary, Abell 2390 shows a larger cooling region, which 
corresponds to: 1. smaller values for /3 and 2. larger for 
L. 

We underline that in our approach it is very important 
to find these kind of correlations: if we want to make the 
scalar field mechanism the most general and basic possi- 
ble, we need to find any possible link among its param- 
eters and the physical properties of our analyzed gravi- 
tational systems, so being able to perform forecasts for 
other scaling properties and to recognize that a scalar 
field is acting even when it is not directly possible to 
derive easy quantities to be compared with data. 

Going on in this search, we have plotted the scalar field 
length L versus the virial radius of any cluster (bottom 



panel in Fig. ([5])), calculated using the relation [l^ IbsI]. 

W = 1.95h-iMpc(^y'' , (74) 

where h = 0.742 [s^] and < T > is the average tem- 
perature for any cluster, derived from gas- weighted fit to 
the total cluster spectrum excluding the central region, 
and resulting more directly correlated to clusters mass 
than the X-ray emission weighted temperatures. Even 
if the two lengths result very different, with the scalar 
field parameter L always smaller than the virial radius, 
we have a phenomcnological relation between them (we 
have performed a weighted linear fit which results better 
than any higher order polynomial), 

log L = 3.0082 • log r„r - 7.0476 . (75) 
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Table II. Clusters of galaxies. Column 1: Cluster name. Column 2: cluster total mass. Column 3: gas mass. Gas and total mass 
values are estimated at r = Vmax- Columns 4: gas mass weighted average temperature. Column 5: virial radius. Column 6: 
minimum observational radius. Column 7: maximum radius. Column 8: coupling parameter /? from scalar field {la confidence 
interval). Columns 9: gravitational length L from scalar field (lu confidence interval). 
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Figure 5. Scalar field length plotted versus mean (gas density weighted) cluster temperature {top panel) and the virial radius 
{bottom panel). Objects in brackets are excluded from fits as described in § (|IVB|I . 



From it one can see that the scalar field length scales as 
L (X r^j^ ~ Myir, so that it may possible to find a relation 
between scalar field parameters and the virial (Navarro - 
Frenk - White (NFW) model based) ones. 
We have also checked a possible relation between L and 
the mean temperature < T > (top panel in Fig. (O), 
with the final result: 

log L = 1.51085 • log < T > +1.79122 . (76) 

This last results is again very intriguing, because for clus- 
ters of galaxies we have a mass-temperature relation [l^ , 

Ma/T3/2 cx Ho/H{z) , (77) 

where A is the overdensity level relative to the critical 



density at the cluster redshift, so that Migo = Afuir- 
Then we have to consider that for most of our clusters 
the redshift is very small so that the ratio Hq/H{z) is 
almost equal to 1. We remind that: 

L(x<T >^/2 and L oc r^^ - M„„ , (78) 

so we can deduce that our scalar field length follows 
perfectly the mass-temperature relation, being L cx< 
T >'^/^cx Mmr and that any scatter in the previous rela- 
tion can be attributed to differences in clusters redshifts. 

Finally, we plotted in Fig. (|6]) the scaled temperature 
profiles versus the distance from the center of any clus- 
ters scaled with respect of the scalar field length L ob- 
tained by fit. It is extremely interesting and not obvious 
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Figure 6. Temperature profiles for all clusters plotted as a function of distance from the center and in units of the scalar field 
length L. The temperatures are scaled to the mean (gas density weighted) cluster temperature. Only the two group galaxies, 
MKW4 and RXJ1159, have not been plotted. Colors in the bottom panel mean: black for clusters with < T > greater than 5 
Kev; cyan for < T > between 2.5 and 5 keV; green for < T > less than 2.5 KeV; grey for Abell 2390. 



that the profiles are all re-scaled and self-similar as it 
happens in the usual approach. We can even say some- 
thing more: in the scalar field approach we can see that 
some properties which are depicted in Fig. (16) of [soj . 
as for example the different profiles among subgroups of 
clusters with mean temperature in different ranges (less 
than 2.5 keV, between 2.5 and 5. keV and greater than 
5 keV) with distances rescaled with respect of r^oo, here 
disappear. In fact, all the clusters form an homogeneous 
gravitational family (the only exception being Abell 2390, 
as discussed before). We could say that the scalar field 
length L contains more fundamental information than 
the virial radius. 



C. LSB galaxies: results 

Taking into account LSB galaxies we have one param- 
eter more than the clusters case: together with the in- 
trinsic scalar field parameters we have the stellar mass- 
to-light ratio, Y,, which wc need for the conversion from 
the stellar surface photometry to the stellar mass den- 
sity. In principle this parameter could be left free; but 
we have decided to put a prior on it derived from litera- 
ture. In [s^ they investigate star formation history and 
chemical evolution of LSB galaxies by modeling their ob- 
served spectro-photometric and chemical properties us- 
ing a galactic chemical and photometric evolution model 
incorporating a detailed metallicity dependent set of stel- 
lar input data. Results show that for this class of 
galaxies has usually values between ^ 0.5 and ^ 2. For 
this reason we have constrained the mass-to-light ratio to 
the interval [0; 5] as a conservative hypothesis, and only 
in one case (UGC3851) we had to enlarge the interval up 
to 10. 

An important assumption about this parameter is that 
we assume it being constant over the whole range of data; 
this is an unavoidable and the most general assumption 



one can do without having details of stellar population 
distribution inside LSB galaxies. 

A comment is in order for a correct evaluation of the re- 
sults: during the analysis we have compared the theoreti- 
cal rotation curves coming out from our scalar field model 
using only the observable-matter densities described in 
[79} with data firstly published in 31 1. These data, avail- 
able in the SIMBAD data base, consist of the contribu- 
tions to the total rotation curve coming from the stellar 
component (with an assumed mass-to- light ratio = 1) 
and the gas one (with gas density normalized with a fac- 
tor 1.4 for taking into account helium contribution), to- 
gether with the total observed rotation curve (which can 
be considered an effect of dark matter in a CDM sce- 
nario or coming out from the interaction with a scalar 
field as in our approach). Both the total and the gas 
rotation curve have been submitted to a smoothing pro- 
cedure for deleting irregularities mainly coming from two 
different elements. First, the assumption that has been 
made when deriving mass models from rotation curves is 
that there is symmetry, that all mass in on circular orbits 
and that there is continuity with radius. Raw data show 
scatter and non-circular motions which can produce vir- 
tual or ambiguous rotation velocities. Second, gas densi- 
ties often show small-scale structures, irregularities and 
look clumpy at any distance from the center, and this can 
give possible large fluctuations in the rotation curve. In 
[3l| the smoothing procedure is tested of course, and the 
smooth rotation curve is a very good approximation to 
the raw rotation curve. But some discrepancies with our 
theoretical estimations can be find because of wc have 
used raw gas density profiles from (7^ . 
If we consider that all gas contributions to rotation curve 
are multiplied to /? only, while the star ones to the com- 
bination /3 • Y*, in some cases, depending on the fit values 
and on their contribution, irregularities in the gas distri- 
bution are emphasized and may affect the total rotation 
curve profile. 
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Table III. LSB galaxies. Column 1: UGC number. Column 2: Distance from the source literature. Column 3: disk central 
surface brightness in R-band, corrected for galactic extinction and inclination. Column 4: disk scale length. Column 5; total 
HI gas mass. Column 6: Maximum rotation velocity. Column 7; best fit stellar mass to light ratio [la confidence interval). 
Column 8: coupling parameter j3 from scalar field (Icr confidence interval). Columns 9: gravitational length L from scalar field 
{la confidence interval). 
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We can now discuss the obtained values for the stel- 
lar mass-to-light ratio. It is easy to see that 10 galaxies 
have values compatible with the prescribed range by j86| ; 
one (UGC3851) has a higher value which is hard to ex- 
plain in term of reasonable population synthesis; and four 
(UGC4325, UGC5721, UGC7603, UGC8286) have values 
smaller than 0.5. Among these, one is particularly prob- 
lematic, i.e. UGC7603, because it has a mass-to-light 
ratio too much small, Y* = 0.017, while the others can 
be again considered acceptable. 

UGC3851 is a real challenge for both our approach and 
the more traditional one [3l| . The very linear rise in the 
inner region rapidly changes in a flat part at larger radii 
so that it is very hard to reproduce this sharp change in 
the slope. There are two possibilities for this behavior: 
the HI curve in outer points is underestimated, or there 
are non-circular motions related to a bar-like structure 
and a star forming region in the center which affect the 
Ha observations. It is interesting to note the many sim- 
ilarities between our analysis and that one in [sij: both 
their NFW and isothermal halo model overestimate ve- 
locity in the inner region, with a larger deviation in the 
NFW model than the isothermal one, and larger in the 
inner region than in the outer one, as in our case. Then 
their derived maximum value for the stellar mass-to- light 
ratio is = 5.4, very close to our value. Anyway, all the 
following relations will be derived without considering 



this galaxy because no definitive conclusion can be done 
about it. 

Among the four galaxies with a low mass to light ratio, 
one of the most problematic case is UGC5721, showing 
a very narrow bi-modality in the parameters distribu- 
tion that cannot be avoided and resolved even chang- 
ing parameters of the trial distribution in the MCMC. 
The convergence test inevitably fails in this case, but we 
have assumed as best fit results the values associated to 
the highest peak in the likelihood distribution shown in 
Fig. ([9]), which produce a very good fit to the observa- 
tional rotation curve. A discrepancy can be found around 
2 kpc but as the large error bar shows, in this region there 
are some observational difficulties related to gas observa- 
tions which produce the detected faster rise in velocity. 

For what concerning the scalar field parameters, the 
coupling constant /3 seems to be very well constrained 
in the range [0.957; 2.416], if we exclude the previously 
discussed case of UGC3851, whose value is very low, 
/? = 0.238, probably because of the degeneracy with Y*, 
and the opposite situation of UGC4325, which exhibits 
a large value, /? = 4.339, and also larger errors on this 
parameter than other galaxies. It is not a causality that 
this is also one of the four galaxies with a stellar mass- 
to-light ratio lower than the simulation limit. 
As in clusters case, the scalar field length shows a larger 
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Figure 7. Rotation curves of LSB galaxies. Dots are velocities from data; solid line is the theoretical model, Eq. (|69|) : dashed 
line is the stellar contribution to the rotation curve assuming F* = 1; dotted line is the gas contribution. 



spread, mainly scaling with the dimension of any LSB 
galaxy. One problematic situation is UGC8837, which 
clearly deviates from the general trend for what concern- 
ing the value of scalar field length L, resulting it very 
small. Even if the fit of our model to data is very good, 
we do not consider anymore this galaxy in the following 
discussion because its results are strongly related to the 
low quality of data. As discussed in |3l|, the Ha are of 
good quality, but do not correspond very well to the HI 
profile, with a large difference in their systemic velocity 
probably due to the interference of inner star formation 
region and there located non-circular motions. A fur- 
ther error source can be also the inclination, being this 
galaxy almost edge-on, the most problematic configura- 
tion in this case. 



Just by visual inspection it is possible to see that 
the general shapes of the rotation curves are quite 
well predicted by our model inside the error confidence 
region, even considering some singular cases where 
are present irregularities coming from gas distribution 
(clumpy peaks, change in the profile convexity). In 
evaluating the goodness of our fitted parameters we 
have to note that all the rotation curves are limited to 
a certain distance from the center of any galaxy. It is 
well known that spiral galaxies show a great variety 
of rotation curves profiles, with flat velocity plateaux, 
different slopes (both increasing and decreasing) in 
outer regions; but in our case the necessity to match 
stellar photometry with more extended gas observations, 
results in a limited distance range which may affect our 
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Figure 8. Same of Fig. (O. 
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Figure 9. Same of Fig. ([Tj for UGC5721. In upper plots are shown tlie marginalized probability distributions for /? and y* 
obtained from MCMCs and the unresolved bi- modality discussed in § pV C[) . 



results. As pointed out in [T^ this could not solve the 
previously described degeneracies betvifeen the fitting 
parameters, or to induce their wrong estimations, so 
requiring more extended data for a more exact analysis. 

Then, for best appreciating the statistical validity of 
our analysis, we can compare them with a recent work 
[sij l on the same class of galaxies but based on MOND as 
alternative theoretical model to dark matter. MOND is 
well-know to be the main and most successful alternative 
scenario in explaining rotation curves of spiral galaxies, 
even if not being satisfactory (when even completely un- 
able) in describing mass profiles in clusters of galaxies. 
In [8l| MOND results to predict quite well LSB rotation 
curves; a look to their figures show that even MOND is 
unable to explain all the features that appear in such pro- 
files, and even some open questions remain, that seem to 
be intrinsic in the model while only weakly depending on 
possible observational sources of uncertainties. 
A positive step of our model, is that we are able to obtain 
sufficiently good fits for two very different gravitational 
systems as clusters of galaxies and LSB and opening at 
the same time the possibility of a unique theoretical back- 
ground underlying both them. 

Moreover, some of the correlations they find between 
MOND acceleration and physical properties of LSB 
galaxies can be found in our approach too. For example, 
wc find a correlation among the parameter /3 and the ex- 



trapolated central disk surface brightness Hq^r, even more 
general than their on, because in our case it well fit all the 
LSB galaxies, without any cut-of dependence (their cor- 
relation strongly depends on galaxies with higher surface 
densities). Two possible fits are: 

/3 = -0.277167 • /io,fl -h 7.78297 (79) 

and 

log (3 = -3.80827 • log ^o,_r + 5.32519 , (80) 

and, as shown in the top panel of Fig. ([TU)) . they are al- 
most indistinguishable. Then, we also have a correlation 
among both scalar field parameters and the maximum ro- 
tation velocity, while in MOND no correlation has been 
found. In the middle and bottom panels of Fig. pH)) we 
have: 

log /3 = 0.354674 • log Vmax - 0.475189 and (81) 

logL = 3.81376 • \ogVmax - 6.14695 . (82) 

Still some questions have to be better studied, because 
for LSB galaxies we could not have performed a detailed 
analysis of scaling or structural parameters as we did with 
clusters of galaxies. 

Anyway, the selected sample has given important de- 
tails and perspectives about the possibility that scalar 
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Figure 10. Correlations among scalar field parameters and structural properties of LSB galaxies. {Top panel.) Coupling 
constant /3 versus the central disk surface brightness. The dotted line is the linear fit; the dashed line is the log-log fit. (Middle 
panel.) Coupling constant /3 versus the maximum rotation velocity. (Bottom panel.) Scalar field length L versus the maximum 
rotation velocity. Objects in brackets are excluded from fits as described in § (|1V[) . 
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field works well even at these scales. But the large dis- 
persion in scalar field parameters show that something 
more accurate should be done. The used sample of LSB 
galaxies is limited by the fact that many rotation curves 
are not smooth, not extended to large radii; and then 
they constitute a limited sample of a more complex and 
extended class of gravitational systems as galaxies. In a 
forthcoming paper we are going to revisit all these ques- 
tions, enlarging the galaxy sample from dwarf and ir- 
regular galaxies, to high surface brightness (HSB) spiral 
galaxies, and to elliptical galaxies. 

D. Unified picture 

Even with all the previous caveats in mind, it seems 
that results are consistent when we compare clusters and 
LSB analysis. We can say more: even if we will only de- 
rive phenomenological and visual relations, whose phys- 
ical meaning has to be studied in more detail, they can 
help us in finding problematic cases and thus verifying a 
posteriori if they are such for intrinsic problems of our 
model or for something different. 

Figs. PT|) - ([15]) show that a possible general trend 
including both clusters and LSB galaxies is feasible. 
In Fig. pT|) a correlation between the scalar field length 
L and the gas density of any object is possible. We have 
considered only gas contribution, because we do not want 
to take into account dark matter and use only visible 
matter: in clusters of galaxies, gas is largely the main 
contribution to the visible mass; and for a self-consistent 
discussion, we have considered an analogous quantity for 
LSB galaxies too, also considering that for calculating 
the stellar mass contribution, we have to use the mass- 
to-luminosity ratio y*, being this a parameter fit. 

For clusters (without the problematic cases: Abell 
2390, MKW4, RXJ1159) we have: 

logL = -1.56605 • log pgas + 8.99653 ; (83) 

for LSB galaxies (without UGC3851, UGC4325, 
UGC8837): 

logL = -1.8617 • log Pgas + 9.95973 ; (84) 

while for the total sample (without the previous excep- 
tions): 

logL = -1.85764- log Pgas + 10.1326 (85) 

We can note that a small difference in the slope between 
the total and the cluster-only sample, while the one from 
spiral is practically equivalent. This may lead us to say 
that a possible general trend is present, and that, likely, 
it can be made more clear when adding further interme- 
diate data, as group of galaxies and elliptical galaxies, or 
smaller scales ones. 

The same, but for the coupling constant /3, is shown 
in Fig. (HH). For clusters (without Abell 2390, MKW4, 



RXJ1159) it is: 

log /3 = 0.0946259 • log pgas - 0.0237813 ; (86) 

for LSB galaxies (without UGC3851, UGC4325, 
UGC8837): 

log /3 = 0.446629 • log pgas - 1 -92872 ; (87) 

and for the total sample (without the previous excep- 
tions): 

log /3 = -0.0290609 • log pgas + 0.444376 (88) 

In this case we see that the total-sample slope is differ- 
ent from the other two cases, even if the one from the 
cluster-only sample is very small, and also by visual in- 
spection clusters seem to be less spread around the main 
general relation than LSB galaxies. These last ones, on 
the contrary, seem to show a proper intrinsic slope, even 
if it can depend on the previously described problems we 
arc facing with when working with a restricted galaxy 
sample, or with not enough extended rotation curves. 

We note that it is much important to verify a very low 
value (« 0) for /3 because we remind that in ? III Al we 
assumed that /3 is constant or at least has a weak depen- 
dence on scale. This hypothesis is partially confirmed by 
clusters and by the total sample fit. 

The same conclusion can be derived from Fig. ([13]) 
where we plot f3 versus L: it is evident that there is quite 
no dependence for (3 on the gravitational scale. In this 
case we also note that the global fit feels more clusters 
than LSB, having performed a weighted fit and clusters 
showing best constraints on scalar field parameters. For 
clusters (without Abell 2390, MKW4, RXJ1159) we have: 

log/3 = -0.11781 • logL + 0.697316 ; (89) 

for LSB galaxies (without UGC3851, UGC4325, 
UGC8837): 

log/3 = -0.00355508 • logL + 0.197588 ; (90) 

and for the total sample (without the previous excep- 
tions): 

log /3 = 0.0176545 ■\ogL + 0.306458 (91) 

When plotting the scalar field parameters versus the 
total gas mass enclosed in the considered gravitational 
structures it is more evident the need of more objects 
for giving more detailed and best constrained results. In 
Figs. ([T4]) - (fTS)) it is evident the big void lying among 
clusters and LSB regions. 

We have given possible fits even in this case; for clus- 
ters (without Abell 2390, MKW4, RXJ1159) we have: 

log L = 0.558047 • log + 0.273985 , (92) 



log /3 = -0.0745605 • log + 0.702031 ; (93) 
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Figure 11. L vs pgas- The dotted lines are the singular fits to clusters and LSB samples. The dashed line is the fit to the 
total (clusters + LSB) sample. The fits are weighted with errors on parameters derived from MCMCs. Objects in brackets are 
excluded from fits as described in § (|IV|) . 
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Figure 12. j3 vs pgas- The dotted lines are the singular fits to clusters and LSB samples. The dashed line is the fit to the 
total (clusters + LSB) sample. The fits are weighted with errors on parameters derived from MCMCs. Objects in brackets are 
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while for LSB galaxies (without UGC3851, UGC4325, 
UGC8837): 



logL = 1.141 • log 



109 Afp 







0.592915 , (94) 



M 

log p = -0.355447 • log ^ + 0.193349 ; (95) 



109 Af( 



o 

For the total sample (without the previous exceptions) 
we have tried two different fits for the length L: 

logi = 1.141 -log ^^^+ 0.592915, (96) 



M 

log L = 3.67182 • log ( 1.38315 log ^ | + 1.59465 

(97) 

but without intermediate data we cannot infer any con- 
clusion; while for /3 we have: 

M 

log /3 = 0.0335265 • log „^°" + 0.198459 . (98) 



V. CONCLUSIONS 

In this work we studied the dynamical properties 
of several astrophysical systems within the theoretical 



framework of scalar theories. We investigate whether 
there are evidences for a scalar field in the considered 
astrophysical systems and if it is possible to observation- 
ally detect it. We have taken into account three different 
classes of objects: supernovae, low surface brightness spi- 
ral galaxies and clusters of galaxies. Results show that: 
i) there is an intrinsic difficulty in extracting information 
about scalar field mechanism (or more generally about a 
varying gravitational constant) from supernovae; ii) a 
scalar field can fairly well reproduce the matter profile 
in clusters of galaxies, estimated by X-ray observations 
and without the need of any additional dark matter; Hi) 
good fits to the rotation curves of low surface brightness 
galaxies, using visible stellar and gas mass components, 
are obtained. 

These results show that different astrophysical system 
can be used as different tracers of the same physical 
mechanism. Moreover, they point towards the possibility 
of a unifying view of dark matter and dark energy via a 
scalar field, at least at galactic and cluster scales [2^ . 
The main criticism of the approach is related to the fact 
that the very different physical properties and evolution 
of the considered astrophysical systems could insert un- 
wanted biases and priors leading to a wrong overall pic- 
ture of the underlying cosmological model. This short- 
coming could be partially avoided if homogeneous and 
well calibrated samples of data at low, medium and high 
rcdshifts will be achieved in future. 
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